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ABSTRACT 

We report on a survey for narrow (FWHM < 600 kms -1 ) Civ absorption lines m a 
sample of bright quasars at redshifts 1.8 < z < 2.25 in the Sloan Digital Sky Survey. 
Our main goal is to understand the relationship of narrow C iv absorbers to quasar 
outflows and, more generally, to quasar environments. We determine velocity zero- 
points using the broad Mgn emission line, and then measure the absorbers' quasar- 
frame velocity distribution. We examine the distribution of lines arising in quasar 
outflows by subtracting model fits to the contributions from cosmologically intervening 
absorbers and absorption due to the quasar host galaxy or cluster environment. We 
find a substantial number 43 ±6 per cent) of absorbers with Wq 1548 > 0.3 A in the 
velocity range +750 kins" 1 < v < +12000 kms -1 are intrinsic to the AGN outflow. 
This 'outflow fraction' peaks near v = +2000 kms -1 with a value of /outflow — 
0.81 ±0.13. At velocities below v w +2000 kms -1 the incidence of outflowing systems 
drops, possibly due to geometric effects or to the over-ionization of gas that is nearer 
the accretion disk. Furthermore, we find that outflow-absorbers are on average broader 
and stronger than cosmologically-intervening systems. Finally, we find that ~ 14 per 
cent of the quasars in our sample exhibit narrow, outflowing C iv absorption with 
T/T/A1548 > 0.3A, slightly larger than that for broad absorption line systems. 

Key words: quasars: absorption lines - galaxies: active - intergalactic medium - 
accretion, accretion discs 



1 INTRODUCTION 

High speed outflows observed in rest-frame UV absorp- 
tion lines are common components of active galactic nu- 
clei (AGN). The ejected gas is believed to be driven off of 
the accretion disk around a central super-massive black hole 
(SMBH) by either radiation pressure (Proga et al. 2000 and 
2004) or magneto-centrifugal 'forces' (Everett 2005). In lu- 
minous quasars, the outflows might be sufficiently energetic 
to provide important feedback to the galactic surroundings, 
potentially disrupting star formation and halting further 
growth of the SMBH (Hamann & Ferland, in prep., Hamann 
et al. 2007, Di Matteo et al. 2005, and refs. therein). Quasar 
outflows might also contribute to the blowout of gas and dust 
from young galaxies, and thereby help to distribute metals 
to the surroundings and reveal the central accreting SMBH 
as an optically visible quasar (Silk & Rees 1998, Moll et al. 
2007). 

Nonetheless, AGN outflows remain poorly understood. 
Most quasar studies have focused on the blueshifted broad 



absorption lines (BALs), which very clearly identify pow- 
erful outflows at typical velocities of 5000 to 20000 kms -1 
(Korista et al. 1993, Weymann et al 1991). BALs are de- 
tected in ~10 per cent of optically-selected quasars (Trump 
et al. 2006), but the flows themselves might be ubiquitous 
if, as expected, the absorbing gas subtends just a fraction 
of the sky as seen from the central continuum source (e.g., 
Hamann et al. 1993). 

There is, however, a wider range of observed outflow 
phenomena than just the BALs. For example, there are the 
so-called 'associated' absorption lines (AALs) that appear in 
roughly 25 per cent of quasar spectra (see also Vestergaard 
2003). These features are much narrower than the BALs, 
with full widths at half minimum (FWHMs) of several hun- 
dred kms -1 or less, and they have redshifts that are, by def- 
inition, within a few thousand kms -1 of the emission line 
redshift, i.e., z a bs ~ z qua3a r (e.g., Foltz et al. 1986, Anderson 
et al. 1987). However, unlike the BALs, the identification of 
AALs with quasar outflows is often ambiguous. Weymann 
et al. (1979) noted that AALs can have a variety of origins, 
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including quasar host galaxies or cluster environments, and 
cosmologically intervening (unrelated) material, in addition 
to quasar outflows. Nonetheless, there is growing evidence 
that a substantial fraction of AALs and other narrow ab- 
sorption lines (NALs) at even larger velocity shifts do form 
in quasar outflows. In their sample of 66 quasar spectra, 
Weymann et al. measured a strong excess of C IV systems 
with Zabs — z qua ,sar compared to z abs « z quasar . Although 
the statistics were somewhat poor, they measured a quasar 
rest- frame velocity width for this excess of « 1200 km s _1 , 
centred near v — 0, which they attributed to the quasar 
cluster environment, and a significant tail extending out to 
blueshifts of ~ 15000 km s _1 which they identified with an 
ejected/outflow component. 

More recently, Richards (2001) found a ~36 per cent ex- 
cess of narrow Civ absorbers with v > 5000 kms -1 in flat- 
spectrum quasars as compared to steep-spectrum quasars. 
This was interpreted as evidence that over a third of high- 
velocity NALs are ejected from quasars. Well-studied in- 
dividual cases provide direct evidence for an outflow ori- 
gin of some NALs based on, for example, variable line 
strengths, resolved line profiles that are smooth and very 
broad compared to the thermal speeds, and doublet strength 
ratios showing that the absorber covers only part of the 
background light source along our observed sightlines (e.g., 
Hamann et al. 1997a and 1997b, Barlow & Sargent 1997, Mi- 
sawa et al. 2007, and refs. therein). Two variability studies, 
in particular, indicate that at least ~21 per cent of quasar 
AALs form in the inner, dynamic regions of quasar outflows 
(Narayanan et al. 2004, Wise et al. 2004). 

Determinations of the detection frequency and the ve- 
locity distribution of ejected NALs are essential to our un- 
derstanding of the basic structure and physics of quasar out- 
flows. These results will also impact studies of intervening 
systems. Cosmologically intervening C IV absorbers are used 
to study galaxy halos and the enrichment history of the IGM 
(e.g., Songaila, 2001; Ryan-Weber, Pettini & Madau, 2006; 
etc.) Determining the contamination of intervening samples 
by intrinsic absorbers is therefore pertinent to those studies. 

With the goal of mapping the incidence of narrow C IV 
absorption in quasar rest-frame velocity space, we make use 
of the large publicly available spectroscopic database from 
the Sloan Digital Sky Survey (SDSS; York et al., 2000). In 
§2 we describe our data set. In §3 we describe the analyses 
of the absorber catalog and possible systematic effects, in 
§4 we discuss the implications of our results and in §5 we 
summarize our conclusions. 



2 DATA SET 

2.1 The Quasar Sample 

We selected quasar spectra from the 4th data release of the 
SDSS. The spectra cover a wavelength range of 3820A (with 
resolution R = 1800) to 9200A (R = 2100), corresponding to 
Civ redshifts of z ab3 > 1.468 and resolutions of 143 kms -1 
^ Sv ^ 167 km s _1 . We wish to study absorbers over a large 
span of quasar-frame velocities, from negative velocities (i.e., 
redshifted with respect to the quasar) to large (^ 71000 
kms -1 ) positive velocity corresponding to the red-edge of 
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Figure 1. The redshift distribution of our quasar sample. 



the Lya emission. 1 To cover this complete range requires 
Zquasar > 2.14. However, in order to accurately measure the 
quasar redshifts (see §2.2 below), we require z quasar < 2.25 
so that the quasar Mg II emission line is covered by the data. 
Furthermore, the number of high signal to noise ratio (S/N) 
spectra available decreases with increasing redshift. There- 
fore, to maximize the sample size while balancing velocity 
coverage against S/N, we adopt a minimum quasar redshift 
of Zquasar — 1-8. We select the 1000 quasars with the high- 
est median S/N in the r'-band that meet the above redshift 
criteria, as returned by the DR4 spectroscopic query form. 
This selection results in a sample with r'-band S/N ^ 17.9 
in each spectrum. 

The incidence of absorbers at relatively low quasar- 
frame velocity is of particular interest. Therefore, in order 
to increase the sample size at low- velocity, we select an ad- 
ditional 500 spectra with 1.6 ^ z qua sar ^ 1-8. These had 
r'-band S/N ^ 18.7 in each spectrum. The emission redshift 
distribution of the total sample of 1500 quasars is shown in 
Fig. 1. 



2.2 Velocity zero-point determination 

As quasar broad emission lines are typically blueshifted from 
the systemic (centre-of-mass) redshift of the AGN, with the 
amount of the shift varying from line to line (e.g., Tytler & 
Fan 1992; Vanden Berk et al. 2001), it is a non-trivial mat- 
ter to determine accurate velocity zero-points in individual 
spectra. The SDSS determines quasar redshifts using an au- 
tomated algorithm (Stoughton et al., 2002) to perform a 
weighted matching of detected emission lines to empirically- 
determined wavelengths from the composite spectrum of 
Vanden Berk et al. These empirical wavelengths, which gen- 
erally differ from laboratory wavelengths, are determined as- 
suming that [O ill] represents the systemic redshift. Accord- 
ing to Schneider et al. (2002), the statistical error (measured 



1 Absorber-quasar velocity differences greater than ~ 71100 
kms -1 place the Civ A1548 absorption feature in the Lya forest, 
making identification and measurement prohibitively difficult. 
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from the scatter from the individual lines or the height and 
width of the cross-correlation function) in this method is 
less than S z ~ 0.01 (5 V ~ 1000 kms -1 ) RMS for non-BAL 
quasars and 8 Z ~ 0.01-0.03 (S v ~ 3000 kms -1 ) for BAL 
quasars. Additionally, there are also very likely to be sys- 
tematic errors that depend on which lines are covered, and 
thus depend on redshift in a non-straightforward manner. 
Thus, because accurate velocity zero-points are important 
to this study, the SDSS-provided redshifts are insufficient 
for our purposes. 

It has been shown, however, that the Mg II broad emis- 
sion line has both a relatively small average offset from 
[Om] ((A-UMgii) ~ ± 100 kms -1 ) and a relatively small 
quasar-to-quasar scatter in the offset (uav < a few hundred 
kms -1 ) compared to other broad emission lines. The use of 
Mg n-derived redshifts should therefore reduce the disper- 
sion between the measured and true values of the quasar- 
frame absorption velocities and thus increase the precision of 
our measurements. Even more important than reducing the 
scatter in the measured velocities, however, is to remove any 
systematic shift by using an accurate and appropriate aver- 
age velocity-offset of the Mg II emission line relative to [O in] . 
Tytler & Fan give a value for this offset of (AuMgii) = +101 
kms - with a dispersion of <ta« = 92 kms - , where positive 
velocity indicates a blueshift with respect to [O in]. Richards 
et al. (2002), however, find (AwMgii) = —97 kms -1 and 
(7a, = 269 kms -1 . The Richards et al. results are derived 
from a large sample and, like the present study, using SDSS 
data, but with quasars at significantly lower redshift than 
used here. Mcintosh et al. (1999) and Scott et al. (2000) 
provide data at higher redshift using near-infrared spectra; 
we made use of the nineteen high- 2 values from these two 
sources and find (AuMgii) = +102 kms -1 , but with a large 
dispersion. Removing one outlier reduces the dispersion and 
gives (A-UMgii) = — 1 kms -1 . 

Clearly there is no consensus value for either the offset 
or the quasar-to-quasar dispersion in the offset. Thus, we 
conduct all of our analyses using four different techniques for 
determining quasar redshifts (and thus velocity zero-points): 
Mg n-determined redshifts using Mgn zero-point offsets of 
kms -1 (i.e., no Mgn velocity-offset from [Om]), +102 
kms -1 and —97 kms -1 , as well as the SDSS redshift val- 
ues. We adopt the value (Ta» = 269 kms -1 from Richards 
et al. as our zero-point uncertainty, which represents the 
true scatter in the Mgll-[0 ill] velocity offset convolved with 
the uncertainty in the ability to fit the Mgn emission in 
the SDSS data. Our interactive fitting (described below) is 
likely superior to the automated routine employed by the 
SDSS pipeline. However, comparison of our values to those 
from the SDSS database indicate the improvement is in most 
cases small, and thus the adopted value of (ta» = 269 km s -1 
is likely a fair estimate of our zero-point uncertainty. 

To compute the Mg 11-determined redshifts, we fitted 
a Gaussian to the Mgn emission line for each quasar in 
our sample, assuming the given velocity zero-point offset. 
The fits made use of only the central peak of the emission 
line in order to avoid effects from line-asymmetries. In prac- 
tice, the extent of the fitted region varied, depending on 
the emission line profile. The precision with which we were 
able to determine the peak wavelength also depended on 
the profile. The formal (i.e., fitting) errors were typically < 
10 kms -1 . While the actual Mgn zero-point uncertainties 
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Figure 2. The velocity zero-point correction distribution of our 
quasar sample for ^A-UMgll) = kms -1 . Positive velocities 
correspond to our Mg Il-corrected redshift being larger than the 
SDSS redshift. The distribution has a mean of +582 kms -1 and 
a standard deviation of 684 kms -1 . 



are necessarily larger than this, they are in most cases al- 
ways well below the dispersion as described above. Such fits 
were possible in 92 per cent of the sample; while the Mgn 
emission in the remaining 8 per cent was too weak and/or 
noisy (e.g., due to poorly subtracted night-sky lines) to de- 
termine an accurate fit. For these systems, we applied an 
assumed zero-point correction equal to the average correc- 
tion of the Mgll-fitted sample. The distribution of velocity 
zero-point differences (compared to the SDSS values) for 
the re-determined redshifts is shown in Fig. 2 for no Mgn 
velocity zero-point offset. Positive velocities correspond to 
our Mgll-derived redshift being larger than the SDSS red- 
shift. The distribution has a median value of +545 kms -1 , 
a mean of +582 kms -1 and a standard deviation of 684 
kms -1 . The velocity zero-point difference distribution for 
the SDSS-provided redshifts and those determined from the 
Mgn emission fits in the SDSS database is almost identical 
to Fig. 2, except for a large spike at +1500 kms -1 as the 
SDSS line-fitting algorithm does not allow shifts larger than 
this value. 



2.3 Identification and measurement of C IV 
absorbers 

We fit pseudo-continua (including the broad emission lines) 
to each spectrum and identify C IV absorber candidates us- 
ing software described in Nestor, Turnshek & Rao (2005), 
updated and adjusted to apply to C iv. Each continuum-fit 
was visually inspected and, when necessary, interactively re- 
fit. Approximately 2 per cent of the quasars in our sample 
were either misclassified stars or exhibited BAL absorption 
strong enough that there was little unabsorbed flux in the 
spectral regions of interest and accurate continuum-fitting 
proved impossible. These objects were removed from the 
sample. 

The software used to identify C IV candidates is very 
sensitive to narrow absorption; therefore we are confident 
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that all narrow absorbers above our significance level thresh- 
old 5cr for the A1548 line and ^ 3<j for A1551, see 
Nestor et al. 2005) were found. Each candidate was visually 
inspected to remove redundant and false detections. Care 
was taken to avoid misclassifying the low-ionization pair C I 
A1277 and Ol A1302, which has a similar separation to the 
Civ doublet, as high-velocity (v > 48000 kms -1 ) Civ ab- 
sorption. Furthermore, the Lya and Nv blended emission 
lines occasionally led to difficulties for our continuum-fitting 
software. Thus, while we did not exclude these regions when 
detecting and measuring absorbers, we conservatively re- 
strict the analyses in this paper to velocities v < 60000 
kms -1 in order to avoid potential biases. 

The large majority of the C IV absorbers that we iden- 
tified in our sample were straightforward to measure. How- 
ever, there are also many systems exhibiting broad, resolved 
profiles. Furthermore, as C IV absorption is relatively com- 
mon in quasar spectra, blending of systems at similar red- 
shifts is occasionally a problem in the medium resolution 
SDSS spectra, particularly at low-velocity (z a bs ~ z quaS ar) 
where there is an excess of absorbers (see below). While 
highly accurate rest equivalent width (Wo) and FWHM mea- 
surements are not crucial for the purposes of this paper, 
proper counting of systems is very important. It is typical 
in studies of intervening absorbers to designate a minimum 
velocity separation for which two absorbers are counted sep- 
arately. However, series of systems at low-velocity can be 
blended across large velocity ranges, and narrow systems 
can be partially blended with broad systems, making the ef- 
fective separation for which it is possible to discern separate 
systems widely variable. Therefore, we took the following ap- 
proach, that will be described in greater detail in Rodriguez 
Hidalgo, Hamann & Nestor (in prep.) To each system we fit 
a Gaussian pair having the correct C IV doublet separation, 
the same FWHM, and a A1548 : A1551 doublet ratio con- 
strained to the range 1 : 1 to 2 : 1. When absorbers appeared 
to have distinct but blended components, we simultaneously 
fit multiple pairs to the entire blended profile, with each 
component fit with a single Gaussian doublet pair and con- 
sidered as a distinct absorption system. In order to decide 
the number of components to fit, we always used the mini- 
mum number justifiable by the data; that is, we only added 
an additional component when there was clear evidence for a 
separate absorption feature (member of a doublet). In some 
cases, broad components were not well described by a single 
Gaussian pair, but the data did not justify deblending into 
multiple features according to our counting rules. Nonethe- 
less, a single Gaussian pair per component always provided a 
reasonable estimate of both Wo and FWHM for the system 
in question. Since the problematic systems are exclusively 
broad, they do not enter into the analysis described in this 
paper. 



3 ANALYSIS 

3.1 Defining the narrow C IV sample 

For the present study, we want to exclude the more ob- 
vious outflows lines, the BALs and mini-BALs (§1). We 
therefore removed from our sample all absorbers with mea- 
sured FWHM > 600 kms -1 . We also removed all absorbers 
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Figure 3. Scattcr-plot of W,^ 1548 versus quasar-frame velocity 
for our narrow (FWHMsC 600 kms -1 ) absorber sample. Veloci- 
ties shown are for the Mg o-corrected redshifts assuming no Mg II 
velocity zero-point offset. 

blended with absorbers having FWHM > 600 kms -1 as 
their identification can often be ambiguous. This helps re- 
duce uncertainties associated with deblending and counting 
but, in the end, has little impact on our results because 
the number of such systems is small (<1 per cent of the 
sample). Finally, several spectra in the sample contain ex- 
tremely strong BAL absorption which removes flux from a 
large swathe of continuum. We computed both the BALnic- 
ity Index (BI; Weymann et al. 1991) and the Absorption 
Index (Al; Hall et al. 2002) to quantify these absorbers (see 
Rodriguez Hidalgo et al. for more discussion). After visual 
inspection, it was decided to remove spectra with AI> 3120 
and/or BI> 1250, which totaled 62 quasars. We note that 
this is not the equivalent of removing BAL quasars; rather, 
we remove only those with little 'usable' unabsorbed contin- 
uum (although, for the BAL quasars kept in the sample, we 
do not include narrow C IV systems that are blended with 
BALs.) 

The final sample includes 1409 quasars, which provide 
a catalog of 2566 narrow Civ absorbers with 0.1 < Wq 1548 
< 2.8 A and -3000 < v aba < 70000 kms" 1 , 2009 of which 
have W A1548 > 0.3A. We note that W A1548 and FWHM 
are correlated, so the removal of FWHM > 600 kms -1 sys- 
tems biases our catalog against very strong systems. Figure 
3 shows a scatter-plot of Wq 1548 versus quasar- frame veloc- 
ity for the final narrow C IV catalog. 

3.2 Computing the completeness correction 

Our data set is complete in neither Wq 1548 or velocity-space. 
It is therefore necessary to compute a two-dimensional com- 
pleteness correction. We follow a procedure similar to that 
described in Nestor et al. (2005), whereby we determine 
the minimum W A154S , W m ;„, detectable at our imposed 
significance-level threshold at every pixel and half-pixel in 
each spectrum. We then construct an array on a grid of ve- 
locity and W A154S steps with values corresponding to the 
number of spectra for which a C IV A1548 line of the given 
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Figure 4. The number of usable spectra as a function of veloc- 
ity corresponding to W Q A1548 values of O.lA, 0.2A, 0.3A, 0.4A, 
0.5A and 0.6A, using Mg Il-corrected redshifts assuming no Mg II 
velocity-offset. The 0.5A and 0.6A curves are nearly identical - 
above WJ* 1548 :^ 0.5A our sample is complete in Wq 11548 at all 
velocities. 



Wq 1548 could have been detected at the given velocity, ac- 
counting for the removed regions of spectra having absorbers 
with FWHM > 600 kms -1 . This 'completeness array' can 
then be used to correct any computed incidence of absorbers 
in both Wq 1548 and velocity space. In Fig. 4 we show the 
number of usable spectra as a function of velocity for dif- 
ferent minimum W A154S values. The 0.5A and 0.6A curves 
are nearly identical - above Wq 1548 ~ 0.5A our sample is 
complete in Wq 1548 at all velocities. Although the mini- 
mum detectable W^ 1548 depends, in principle, on the ab- 
sorption FWHM, in practice it is computed assuming unre- 
solved lines. However, any detection bias arising from this 
would only be relevant for lines with Wq 1548 ~ Wmin that are 
also resolved in velocity. However, lines weaker than « 0.5A 
(the strength above which our survey is complete) are gen- 
erally unresolved in the SDSS spectra. 



3.3 Results 

3.3.1 The quasar-frame velocity distribution 

Using the completeness array together with the absorber 
catalog, we can compute the incidence of absorbers as a 
function of velocity, where (3 — v/c. The magnitude of 
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-gp represents the average number of absorbers per kms - 
times the speed of light. Thus the average number of ab- 
sorbers per spectrum can be found by integrating ^5 over 
the relevant velocity range: 



(n) = - 
c 



We show % 

9/3 



dn , 
for W A 



(1) 



^ 0.3A, using Mg il-corrected 
redshifts and assuming no Mgn velocity-offset, in Fig. 5. 
A large excess at v « (i.e., z a b s — Zquasar) is clearly 
detected. The observed distribution includes contributions 
from: (i) cosmologically-'intervening' systems; (ii) 'environ- 
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Figure 5. The incidence of C IV absorbers versus quasar-frame 
velocity for systems with Wq 1548 ^ 0.3A, using Mg Il-corrected 
redshifts and assuming no Mgn velocity-offset. Horizontal bars 
indicate the velocity bins and vertical bars the 1 a uncertainty 
from counting statistics. The curves are maximum-likelihood fits 
to the data in the intervals v < and 40000 kms -1 < v < 
60000 kms -1 , representing the 'environmental' and 'intervening' 
components, respectively (see text). 



mental' absorption which arises in either the ISM of the 
AGN host galaxy or the IGM of the galaxy's group/cluster; 
and (iii) 'outflow' systems that are ejected from the central 
AGN (e.g., accretion-disk outflows). We model the inter- 
vening contribution to the distribution with a step function 
that breaks at v = 0, convolved with the redshift uncer- 
tainty (modeled as a Gaussian with a z = 269 kms -1 ). We 
model the environmental contribution with a Gaussian cen- 
tred at v — 0. To obtain the best-fit parameters, we employ a 
maximum likelihood fit to the unbinned data (absorber cat- 
alog plus completeness array), using only data with v < 
and 40000 kms -1 < v < 60000 kms -1 to constrain the 
fit, where the outflow fraction is assumed to be insignificant 
at these velocities. The physical velocity dispersion of the 
AGN-host's local environment can be estimated from the 
environmental-component by deconvolving the Gaussian fit 
with the assumed redshift uncertainty (269 kms -1 ) and, fol- 
lowing Weymann et al., divide by y/2 to account for the fact 
that what we measure is the difference between the velocity 
of the AGN-host within the environment and that of the 
other halos (although the appropriateness of this last step 
is limited by the presence of absorption in the AGN host 
galaxy, which is not virialized in the local environment). 

The parameters describing the fits are shown in Table 
1 for each zero-point determination method (§2.2). The first 
column indicates the velocity zero-point offset (for the Mg n- 
determined redshifts cases), and the parameters are defined 
by: 



= Nenv X exp(— v /2a env ) 7 



d/3 

V2 x a' env = \fal n 
and 
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Table 1. Environmental and intervening fit-parameters. 



(AuMgii) 
(kms" 1 ) 


-^V env 


(T env 

(km s" 1 ) 


ff enu 

(kms -1 ) 


inter 





31.5 


681 


442 


6.4 


+102 


26.8 


665 


430 


6.4 


-97 


35.4 


694 


452 


6.4 


SDSS 


25.9 


1319 


913 


6.5 



The a en v values are the Gaussian widths before and after de- 
convolution and the N values are the normalization of the fits 
— i.e., the incidence for intervening systems and the incidence of 
environmental systems at v = 0. 



dn int 



dp 



= Ni, 



(4) 



Our assumption that none of the absorbers at v > 40000 
kms -1 or ii < arise in outflows means that our derived out- 
flow fractions (see below) are strictly speaking lower limits 
at all velocities. For example, Richards et al. claim that as 
many as 36 per cent of C IV absorbers at 5000 < v < 55000 
kms -1 form in outflows, while Misawa et al. (2007) claim 
that the fraction of outflow lines at these velocities is ~ 
10-17 per cent. However, it is unclear what the contribu- 
tion of very high-w systems (i.e., v > 40000 kms -1 ) to the 
Richards result is, and the Misawa et al. sample contains no 
systems at v > 40000 kms" 1 that are classified as a reliable 
narrow-line outflow candidates. Thus, while we cannot di- 
rectly test our high velocity assumption, it is likely that we 
are only slightly, if at all, overestimating the contribution of 
intervening systems. 

The assumption that all v < absorbers are environ- 
mental is potentially less reliable, however. It will be an 
appropriate assumption if there are no 'in-falling' systems 
intrinsic to the AGN flow and all of the individual quasar 
velocity zero-points are accurate (or that there are no out- 
flowing systems with velocities that are small relative to the 
velocity zero-point dispersion). However, our zero-points are 
inaccurate with an (assumed) RMS of ~ 270 kms" 1 , and 
the systematic velocity-offset of Mgn from the true veloc- 
ity zero-point is not well constrained (see §2.2). Thus it is 
possible that our environmental fits are also accounting for 
some low-velocity outflow systems. We investigate this issue 
below. 

The uncertainty in the appropriate velocity zero-point 
has a direct effect on the measured velocity-space distribu- 
tion of absorbers. In order to determine the magnitude of 
this potential bias and investigate its affect on the computed 
outflow fraction, we computed §| as well as the maximum- 
likelihood fits using all four of the quasar redshift determina- 
tion methods discussed in §2.2. Fig. 6 shows the low-velocity 
region of ^ computed using: Mg n-determined redshifts 
with (a) no Mgn velocity zero-point offset (i.e., (AuMgii) = 
0), (b) (AwMgii) = +102 kms" 1 , (c) (AuMgii) = -97 
kms" 1 , as well as (d) SDSS redshifts. The curves repre- 
sent the data well in the fit-regions, for all of the offsets 
considered. However, the distribution determined assuming 
(AtiMgii) = —97 kms" 1 overpredicts the incidence at v > 0. 
We note that, if the velocity zero-point is (on average) prop- 
erly determined, the incidence for v > should be larger 
than that for v < 0, since the v > data contains, in prin- 
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Figure 6. Same as Fig. 5, employing Mg Il-determined redshifts 
with (a) no Mgn velocity zero-point offset, (b) ^AiiMgii) = +102 

kms" 1 , (c) ^Aii M gll) = -97 kms" 1 , as well as (d) the SDSS- 
provided redshifts. The width of the environmental-component 
determined from SDSS redshifts is unphysically large, and that 
determined assuming /AvMgll) = ~~ 97 kms -1 gives a poor fit 
to the v ~ data. 



ciple, environmental, intervening, and outflowing Civ sys- 
tems, and we assume the environmental-component is rea- 
sonably symmetric about the true velocity zero-point. We 
obtain a value of a' env = 913 kms" 1 using the SDSS red- 
shifts, which is unphysically large unless the quasar environ- 
ments are predominantly in large clusters (e.g., Becker et al. 
2007). However, Outram et al. (2003) show that the cluster- 
ing amplitude of quasars at z ~ 1.4 is similar to that of 
present-day galaxies. The Mg il-determined redshifts lead to 
values of 430 < o' mv < 450 kms -1 for the dispersion, which 
is comparable to the usual division between large groups and 
poor clusters (e.g., Mulchaey, 2000). 

Considering the quality of the fits, we conclude that, 
while we do not explicitly reject any of the Mg n-corrected 
redshift choices, an average Mgn velocity zero-point offset 
of ~ to +100 kms -1 is most appropriate for our sample. 
Except when explicitly stated, we will use (AtiMgii) = for 
the subsequent analyses. We also note that we computed g§ 
for ranges of W / A154S , but found no significant trends. 
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Figure 7. Top: The excess of the data over the sum of the two 
fits, divided by the data, which represents the minimum outflow 
fraction. Bottom: The data minus our fits to the environmental 
and intervening systems, representing ^ for the ejected (i.e., 
outflow) component only. 



3. 3. 2 The outflow fraction and incidence of outflowing 
absorbers 

The top panel of Fig. 7 shows the data minus the sum 
of the two fits, divided by the data, which represents the 
minimum outflow fraction, f outflow (which is considered a 
minimum for the reasons described above). We find that 
f out flow increases strongly from v = to v ~ +2000 kms -1 , 
where it peaks with f outflow ~ 0.81 ± 0.13, and then de- 
creases slowly out to v ~ +12000 kms -1 . Over the range 
where there is significant evidence for a non-zero outflow 
fraction, v ~ +750 kms -1 to v « +12000 kms -1 , we find 
{f out flow) = 0.43+0.06. Over narrower ranges near the peak, 
we find (f outflow) = 0.57 + 0.10 for +1250 kms" 1 < v < 
+6750 kins" 1 , and f outflow ^ 0.72 ± 0.07 for v ~ +1250 
kms -1 totiRi +3000 kms -1 . The outflow fraction decreases 
below v « +2000 and disappears below v w +750, appar- 
ently indicating an effective minimum projected ejection ve- 
locity for narrow C IV systems. It is also possible that sys- 
tematics such as a strong proximity effect for the intervening 
absorbers or limitations in our ability to properly count sys- 
tems contribute to this decrease. We discuss these and other 
possible biases in §3.4. 

In the bottom panel of Fig. 7 we show the data after 
subtraction of the sum of the fits (using (A«M g n) = 0), 
which represents ^ for the outflow-component only. Using 
(AwMgii) = —97 kms -1 (i.e., Mgn emission redshifted from 
the quasar rest-frame) or SDSS-determined redshifts leaves, 
respectively, marginally- and strongly-statistically signifi- 
cant negative, and thus unphysical, residual incidence for the 
outflow-component at some velocities. Using (AvMgii) = 
(i.e., Mgn emission unshifted from the quasar rest-frame), 
the incidence of systems is consistent with zero at low pos- 
itive velocities: in this case, our assumptions that no ab- 
sorbers arising in outflows have measured velocities with 
ii < is likely valid. Attempts to fit the entire span of ve- 
locities (i.e., including < v < 40000 kms -1 , see §3.3.1) 
with, e.g., a power-law or exponential for v > 0, failed since, 



as is clear from the bottom panel of Fig. 7, these were not 
appropriate descriptions of the outflowing component. Con- 
sidering the results employing Mg n-corrected quasar red- 
shifts, we find that the incidence of narrow C IV absorbers 
with Wq 1548 > 0.3 A presumably forming in accretion-disk 
outflows peaks at v ~ 2000 kms -1 and exhibits a wing that 
extends out to at least v ~ 9000 kms -1 , for all (AuMgii) 
values considered. 

3.4 Possible systematics 

In this section we discuss potential biases introduced by 
the method we employ to 'count' systems as well as the 
assumptions behind the modeling of the distribution of 
cosmologically- intervening absorbers. 

As narrow C IV absorbers cluster at low (< 1000 km s -1 ) 
quasar-frame velocities (Nestor, Hamann & Rodriguez Hi- 
dalgo, 2007), removing regions with broad absorption is 
more likely to also remove narrow absorbers, compared to 
randomly selected regions at similar velocity. Since broad 
absorbers are more common at low positive velocity, it is 
a concern that this bias may affect the measured distribu- 
tion of narrow systems. Furthermore, this region of velocity 
space exhibits (in principle) all three of our categories of ab- 
sorbers, thus increasing the chances for blending of outflow 
and non-outflow systems to occur. Therefore, the apparent 
dearth of narrow outflow-component systems at v < 2000 
kms -1 (Fig. 7) could conceivably be in part a byproduct of 
our removal of broad systems and absorbers blended with 
broad systems. To test this, we re-computed |^ without 
rejecting any absorbers based on FWHM or blending with 
broad systems. However, the only qualitative difference in 
our results was an extension of the high-velocity tail of out- 
flowing systems out to v « 25000 kms -1 , indicating that 
broad (FWHM> 600 kms -1 ) absorbers have a larger ex- 
tent in ejection-velocity than do narrow systems. This will 
be discussed in further detail in Rodriguez Hidalgo, Hamann 
& Nestor (in prep.) 

Alternatively, over-zealous de-blending could, conceiv- 
ably, influence both the shape and mag nitude of f^. To 
investigate this concern, we combined all absorbers with ve- 
locity differences less than the sum of the two HWHM val- 
ues plus 150 kms -1 (i.e., the SDSS resolution) into a single 
absorber to be counted once. We then recomputed |^ and 
the fits with this reduced catalog. This resulted in a slight 
decrease in the outflow fraction for bins where it was non- 
zero, and virtually no change in the shape or range over 
which f out flow > is significant. Even when we took the 
extreme step of combining all absorbers within 500 kms -1 
plus the sum of the two HWHM values, the shape of ^ for 
outflow-only systems remained the same, although the mag- 
nitude was reduced by ~ 20 - 25 per cent. Thus, we believe 
that any inconsistencies in our counting method is at a level 
small enough to be safely ignored. 

Another concern is that cosmologically-intervening sys- 
tems may not be uniformly distributed in quasar-frame ve- 
locity for our quasar sample. For example, it is known that 
at high-z the incidence of intervening C IV absorbers de- 
creases with increasing redshift (Sargent, Boksenberg & Stei- 
del, 1988; Misawa et al., 2002). Larger redshift correlates 
with smaller velocity in Figs. 5-7. According to Monier et 
al. (2006), the incidence of W A1S48 -> 0.3A systems show 
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little evolution over the range 1.46 < z < 2.25. Misawa et 
al., however, claim a ~ 50 per cent decrease over the same 
range. While the decrease in f out ji ow is too abrupt to be en- 
tirely caused by redshift effects, we nonetheless investigated 
the maximum possible magnitude of this effect by running 
50 Monte Carlo simulations of intervening systems in our 
data. We randomly distributed absorbers in the spectra of 
our sample in redshift-space using the parameterization 
from Misawa et al., and converted the resulting distributions 
into velocity-space to determine the modelled fe. This re- 
sulted in < 10 per cent decrease in |^ from v « 45000 
kms" 1 to v — 0, or a difference of < 0.7 in fe. 

Potentially more important is the possibility of a strong 
proximity effect, similar to that seen for Lyman-a forest 
lines, causing a decrease in ^ for the intervening systems. 
While the magnitude of any proximity effect likely depends 
on the line-strengths being considered, C IV has a much 
higher ionization energy than does Hi (64.5 eV compared to 
13.6 eV) and thus any C IV proximity effect should be signif- 
icantly weaker than that for H I. Thus, we ran 50 additional 
Monte Carlo simulations of intervening systems using the 
proximity effect results for the Lyman-a forest from Scott et 
al. (2000) , considering this a strong upper-limit to any C IV 
proximity effect. Doing so, we find a linear decrease in ^ 
from no difference at v — 2000 km s _1 to a 50 per cent deficit 
of systems at v = 0. This corresponds to a difference in |^ 
at v = of ~ 3.5. Thus, a C IV proximity effect as strong as 
what is claimed for Lyman-a would only marginally change 
the results described above, and then only at v ~ 0. 

Finally, we also considered the effect of an underestima- 
tion of the velocity zero-point dispersion. A larger dispersion 
would cause more low-velocity intervening systems to scat- 
ter to v < 0. We repeated the modeling of the environmental 
and intervening components doubling the velocity zero-point 
dispersion and found no qualitative changes to our results. 

We are limited in our ability to determine true individ- 
ual velocity zero-points, to deblend ('count') absorbers into 
physically-distinct systems with complete accuracy, and to 
know the true distributions of the non-outflowing absorber 
populations. These systematics almost certainly affect our 
results to some degree. Nonetheless, the maximum estimates 
of the effects from all of the systematics that we considered 
are relatively small and have no significant impact on our 
qualitative results. 



4 DISCUSSION 

4.1 Properties of the outflowing absorbers 

We find a significant detection of narrow C IV absorption 
over that which is expected from cosmologically-intervening 
and host galaxy-environment contributions, putatively oc- 
curring in accretion disk winds, in the quasar-frame velocity 
range 750 kms -1 < v < 12000 kms -1 . Using eqn. 1 to- 
gether with the §| values for the outflow-component (Fig. 
7), we can compute the average number of outflowing ab- 
sorbers per quasar. We find (n utfiow) = 0.25 ± 0.02 using 
(A«Mgii) = 0. The value depends only slightly on the choice 
of (AvMgii) and the limits of integration. 

While it is likely that some W 7 ^ 1548 ^ 0.3k low-velocity 
(< 750 km s~ ) narrow C IV absorbers arise in outflows, the 
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Figure 8. The normalized (left) and cumulative (right) FWHM 
distributions for high- and low- velocity absorber sub-samples. 
The low-velocity sample is clearly broader. Note that the SDSS 
resolution (~ 150 kms -1 ) sets the lower-limit seen in these data. 
Also shown (solid histograms) are the weighted difference distri- 
butions (using the measured outflow fraction, see text), which 
represents the distributions for outflowing-systems only. 



various possible systematic effects discussed above are un- 
able to completely, or even mostly, explain the turnover in 
|^ for outflow-only absorbers below v < 2000 kms -1 - 
it appears rather to be an indication of a minimum pro- 
jected line-of-sight ejection velocity for narrow C IV absorp- 
tion systems in AGN outflows. This empirical phenomenon 
can help to constrain quasar outflow models. For example, 
the deficit of Civ outflow-systems at v < 2000 kms -1 might 
result from a tendency to 'over-ionize' the low-velocity gas, 
e.g., if that gas is typically nearer the ionizing continuum 
source than the high speed absorbers. This can occur in any 
model of the acceleration (radiation pressure or magneto- 
centrifugal forces), but over-ionization at low velocities has 
a direct physical explanation if there is radiative driving. 
In particular, highly ionized gas is more transparent and 
thus more difficult to push via radiation pressure. High flow- 
speeds are attained only in less ionized gas (that might be 
farther downstream in the flow and partially shielded from 
the intense ionizing radiation, e.g., Hamann & Ferland, in 
prep). 

We can also infer the distribution of other outflow- 
absorber properties by comparing the high-velocity (in- 
tervening) absorbers to those at lower velocities where 
fout flow > but the incidence of environmental systems 
is negligible. To optimize such a comparison, we require a 
range of velocity large enough to include a significant num- 
ber of outflowing absorbers, while restricting ourselves to 
velocities where the outflow fraction is high. Fig. 8 shows 
the normalized (left) and cumulative (right) distributions of 
FWHMs for absorbers with 1800 kms _1 <u< 4400 kms" 1 , 
over which the effective outflow fraction is (/ outflow) = 
0.61 ±0.07 and the environmental contribution is small, and 
30000 kms" 1 < v < 60000 kms" 1 , where we assume the 
absorbers to all be cosmologically intervening systems. The 
low- velocity FWHM distribution is clearly broader than the 
high-velocity distribution. A K-S test gives a zero per cent 
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Table 2. Fraction of quasars containing narrow Civ absorbers. 
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Figure 9. Incidence of absorbers as a function of W A1548 , for 
high- and low-velocity samples. Also shown (solid histogram) is 
the difference between the two distributions which represents the 
distribution for outflow-systems only. 



probability that the two sets of data are drawn from the 
same parent distribution. A broader distribution is expected 
if the low- velocity sample contains outflowing absorbers. We 
also show in Fig. 8 the weighted (assuming f outflow = 0.61), 
renormalized difference between the two distributions, which 
should represent the pure outflow distribution and has a 
mean (FWHM) = 337 kms -1 - though we note that the 
SDSS resolution imposes an artificial lower-limit to the mea- 
sured FWHM of » 150 kms -1 . We also compared the distri- 
bution of A1548 : A1551 doublet-ratios (DR) for the two ve- 
locity sub-samples. The doublet ratio distributions are sim- 
ilar except that the low-velocity sample exhibits relatively 
more completely saturated absorbers: 26 ± 3 per cent with 
DR > 0.9 compared to 15 ± 2 per cent at high-velocity, 
which implies that a third (33 ±6 per cent) of outflowing ab- 
sorbers in this velocity range have DR > 0.9. This result im- 
plies that outflow systems contain fewer weak kinematically- 
distinct components that are unresolved in the SDSS spectra 
than do intervening absorbers. Although individual outflow- 
ing narrow C IV absorbers appear similar to non-intrinsic 
systems in low/medium-resolution data, they are on aver- 
age broader and more often completely saturated. 

Similarly, we can construct the Wq 154S distribution for 
the two velocity ranges. This is shown in Fig. 9. The distri- 
butions are normalized such that they represent the differen- 
tial number of absorbers per A per 3 x 10 5 kms" 1 averaged 
over the velocity range being considered. The distribution 
for outflow-only systems is therefore obtained simply from 
the difference of these two distributions and is shown as the 
solid histogram. It is a much flatter distribution, with a very 
large excess of strong systems compared to the intervening 
absorbers. Notably, almost all of the low-velocity narrow ab- 
sorbers with W A1548 > 1.4A can be attributed to outflowing 
systems. These differences are consistent with the broader 
FWHM distribution of outflow-systems in Fig. 8. 
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0.08 
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0.31 
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The values represent the fraction of quasars that exhibit at least 
one (1+) and at least two (2+) narrow W^ 1548 ^ 0.3A Civ ab- 
sorbers in various velocity ranges, as well as the average number 
of absorbers per spectrum. Values for 'all' absorbers are taken di- 
rectly from the data, while the fractions considering only 'outflow' 
absorbers use our fits to account for environmental and interven- 
ing systems, assuming Poissonian statistics. 



4.2 How many quasars exhibit outflowing 
absorbers? 

Although |p allows one to calculate the average number of 
absorbers that will be found in a given range of velocity 
space per quasar spectrum (i.e., Eqn. 1), it does not contain 
information on the fraction of quasars which exhibit narrow 
C IV absorption. We show in Table 2 the fraction of spectra 
that exhibit one or more and two or more narrow C IV ab- 
sorbers with Wq 154S > 0.3A for different velocity ranges, as 
well as the average total number of systems per spectrum, 
(n). 

As we are primarily interested in the fraction of quasars 
exhibiting absorbers that arise in outflows, we use our mod- 
els for the incidence of environmental and intervening sys- 
tems as a function of velocity (§3.3.1 and Table 1) to predict 
the fraction of quasars having a given number of non-outflow 
absorbers for a given velocity range, assuming Poissonian 
statistics. While it is likely that intervening absorbers ex- 
hibit some degree of clustering, any such possible signal is 
at levels undetectable in our data (Nestor et al. 2007) and 
thus the Poissonian assumption is appropriate. The predic- 
tions for non-outflow absorbers can then be disentangled 
from the observed fractions to calculate the fractions for 
outflowing systems only. These fractions are also shown in 
Table 2. In this manner we estimate that approximately 
fourteen per cent of the quasars in our sample exhibit out- 
flows that are detected in narrow C IV absorption over the 
range < v < 12000 kms , where the intrinsic fraction is 
measured to be non-zero. 

It is important to stress that the outflow-only fractions 
are minima, since we may be counting some outflow-systems 
as non-outflow ones, and they do not account for systems 
that we have eliminated due to blending with broad (FWHM 
> 600 kms" 1 ) absorption (§3.1). However, the total frac- 
tions (i.e., those including the non-outflow absorbers) repre- 
sent firm upper limits for the outflow-only fractions. Thus, 
the overall and calculated outflow-only fractions represent 
windows in which the true outflow-only fractions lie. For the 
reasons described in §3.4, we expect that the true outflow- 
only fractions are near the calculated outflow-only values 
presented in this table. 
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5 SUMMARY AND IMPLICATIONS 

Our main new results are the following: 

1) The distribution of narrow (FWHM < 600 kins -1 ) 
Civ absorbers with W^ 1548 > 0.3A is well described by 
three components: a Gaussian 'environmental' component 
centred at v = kms -1 with a deconvolved dispersion of 
« 450 kms -1 , a flat cosmologically-intervening component 
at v > kms -1 with ~ 6.4, and an outflow component 
that extends from v ~ 750 kms -1 to v ~ 12000 kms -1 
(Figs. 5-7). 

2) The outflow-systems account for a large fraction (> 
61 per cent) of all narrow-line absorbers in the velocity range 
1800 kms -1 < v < 4400 kms" 1 . The outflow fraction peaks 
near v ~ 2000 kms -1 at /outflow — 0.8. The true outflow 
fractions might be even higher, and the tail of the distri- 
bution might reach beyond 12,000 kms -1 , if some of the 
high- velocity narrow systems (30000 kms -1 < v < 60000 
kms -1 ) also form in outflows. This possibility is suggested 
by results from previous studies (e.g., Richards et al. 2001; 
Misawa et al. 2007). 

3) The fraction of quasars with NAL outflows between 
and 12,000 kms -1 is >14 per cent. The fraction of quasars 
with two or more NAL outflow systems at these velocities is 
>6 per cent. 

4) Among the narrow-line systems we consider (FWHM 

< 600 kms -1 ), the outflowing absorbers (e.g., at 1800 
kms -1 < v < 4400 kms -1 ) tend to be broader than the in- 
tervening ones (at v > 30000 kms -1 ). The outflow-systems 
also tend to be stronger and they are more frequently satu- 
rated, as measured by W A154S and the C IV doublet ratio. 

5) High resolution spectra are needed to examine the 
flow kinematics, but the measured line widths <600 kms -1 
imply that the projected flow speeds in these NAL outflows 
are characteristically ~10 times larger than the line-of-sight 
velocity dispersions, i.e., u/FWHM ~ 10. 

6) There appears to be a deficit of outflowing NALs at 
low velocities (v < 1500 kms -1 compared to 1500km s -1 

< v < 8000 kms -1 ). 

These results provide important new constraints on 
models of quasar outflows. For example, the fraction of 
quasars with NAL outflows (>14 per cent) is at least 
as large as the fraction of quasars with classic (strong 
and broad) BALs (~10 per cent), as measured in similar 
optically-selected quasar samples (Trump et al. 2006 and 
refs. therein). Moreover, most of the NAL systems in our 
study appear in quasars without BALs, and we specifically 
exclude NALs that are blended with BALs. Therefore, the 
NAL outflows we measure have no obvious relationship to 
the BALs. They might represent a wholly different type of 
outflow phenomenon - one that is perhaps less massive but 
just as common as the BAL outflows. Perhaps there is an 
evolutionary relationship where massive BAL outflows tend 
to dissipate over time to become NAL outflows. Alterna- 
tively, the appearance of NAL outflow lines could be sim- 
ply a different manifestation of a single quasar outflow phe- 
nomenon, where our viewing perspective through the flow 
determines whether we see NALs or BALs (e.g., Elvis et al. 
2000, Ganguly et al. 2001). Perhaps the reality in quasars is 
a mixture of these two scenarios (Hamann & Sabra 2004). 

The detection frequencies of the outflow lines provide 
important constraints on the flow structures. In particular, 



the true fraction of quasars with BAL outflows, /s, is re- 
lated to the detection frequency by /s > 0.10, while the 
true fraction with NAL outflows is /at > 0.14. If there are 
no orientation biases that affect the detection rates, then the 
covering factors of the outflows (as seen from the central con- 
tinuum source) are >0.14//jv for the NALs and ~0.10//s 
for classic BALs. If all quasars have both NAL and BAL 
outflows as part of a single outflow phenomenon, then the 
covering factor of outflow gas is >0.24. These factors become 
even larger if we include systems not counted here, e.g., the 
broader mini-BALs (Rodriguez Hidalgo et al., in prep) and 
the high-velocity NALs (Misawa et al. 2007, Richards et 
al. 2001). Moreover, all of these results based on detection 
frequencies are lower limits if quasars are characteristically 
fainter when viewed along outflow sightlines (Goodrich et 
al. 1997, Schmidt & Hines 1999), or if quasars appear ob- 
scured and are therefore missed by flux-limited optical sur- 
veys when viewed along any range of sightlines. On the other 
hand, if quasars are preferentially obscured along certain 
viewing angles by, e.g., a dusty torus, the above covering 
factors refer only to the unobscured solid-angle. 

Further insight into the flow structures comes from the 
fact that roughly half of the quasars with NAL outflows (6 
of the total 14 per cent) have two or more of these outflow 
features. Thus, the line-of-sight structure is often complex, 
with either multiple streams or multiple 'clumps' contribut- 
ing to the C IV absorption. In addition, the NAL flow speeds 
are in the low-to-medium range of BALs, suggesting a sim- 
ilar point of origin for these outflows, but the large ratios 
of flow speed to line-of-sight velocity dispersion in NALs, 
u/FWHM ~ 10, contrast markedly with the BALs where 
typically u/FWHM ~ 1. It is tempting to identify each NAL 
system with a clump of gas in the flow, but there is no com- 
pelling model for these features. Clumps or filaments formed 
by instabilities at the ragged edge of a BAL flow should be 
short lived. The measured FWHMs of the NALs, when they 
are resolved in the SDSS data at several hundred kms -1 , 
suggest highly supersonic expansion. 

The possible deficit of low- velocity outflow NALs is an- 
other interesting clue, but it also has no straightforward in- 
terpretation. This deficit might imply that the inner flow 
is too ionized for C IV, or that the clumps or streams pro- 
ducing C IV NALs do not form in the inner flow where the 
speeds are low. There might also be orientation effects that 
preclude our detection of these clumps or streams of C IV 
gas near the base of the flow. 

Although more work is certainly needed to resolve these 
questions, NALs represent another important component of 
quasar outflows. Many studies have focused on the BALs be- 
cause they are the most obvious and most easily measured 
of the UV outflow lines. The overall outflow phenomenon 
is clearly much more diverse, however. Even the NALs and 
BALs together do not tell the whole story (Hamann & Sabra 
2004). These features are at opposite extremes in terms of 
their FWHMs, but in between is a significant population 
of 'mini-BALs' which have intermediate FWHMs and sim- 
ilar (or perhaps characteristically higher) outflow velocities 
compared to the NALs and BALs (Rodriguez Hidalgo et al., 
in prep). The challenge now is to achieve a better physical 
understanding of this full range of quasar outflow features. 
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